Dynamics of Turing patterns under spatio-temporal forcing 
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We study, both theoretically and experimentally, the dynamical response of Turing patterns to 
a spatio-temporal forcing in the form of a travelling wave modulation of a control parameter. We 
show that from strictly spatial resonance, it is possible to induce new, generic dynamical behaviors, 
including temporally-modulated travelling waves and localized travelling soliton-like solutions. The 
latter make contact with the soliton solutions of P. Coullet Phys. Rev. Lett. 56, 724 (1986) 
and provide a general framework which includes them. The stability diagram for the different 
propagating modes in the Lengyel-Epstein model is determined numerically. Direct observations 
of the predicted solutions in experiments carried out with light modulations in the photosensitive 
CDIMA reaction are also reported. 
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The study of pattern dynamics under external forc- 
ing provides a powerful tool to deeply probe their inher- 
ently nonlinear mechanisms under non-equilibrium con- 
ditions. A great deal of attention has been focused on 
resonances or locking of spatially structured states, ei- 
ther stationary or oscillatory, under temporal (spatially 
uniform) OLE, El EL El or spatial (steady) modulations 
[il B M, 111 LLil • Steady patterns in reaction-diffusion sys- 
tems t ypically arise from the celebrated Turing mech- 
anism |f]J, |l2(. According to it, inhomogeneous distri- 
butions of chemical concentrations self-organize sponta- 
neously out of a non-structured medium as a result of 
a competition between autocatalytic reaction steps and 
the differential diffusivities of an activator (smaller) and 
an inhibitor-like (larger diffusion) species. Turing pat- 
terns are endowed with an intrinsic wavelength, depend- 
ing only on the kinetic and diffusion parameters, but lack 
an intrinsic frequency, in contrast to oscillatory chemi- 
cal systems Genuine Turing patterns were first ob- 
served in quasi two-dimensional gel reactors (pre-loaded 
with appropriate chemical indicators) in the CfMA 0] 
and the CDIMA 01 reactions, and appear as patterned 
distributions of iodide. The CDIMA reaction has the 
interesting feature of being photosensitive 0] . 

It seems thus timely to search for generic behavior in 
the unexplored perspective of spatio-temporal forcing of 
pattern forming systems. Specifically, we aim at studying 
the dynamical response of photosensitive Turing patterns 
to the simplest external spatio-temporal forcing consist- 
ing of a travelling- wave modulation of the control param- 
eter associated to the illumination. Through the mecha- 
nism of pure spatial resonance an external frequency will 
thus be imposed in an otherwise nonoscillatory system. 
As a consequence, new nontrivial dynamical modes are 
expected to arise which allow to connect the two trivial 
limiting cases, namely: a travelling pattern locked to the 
forcing at low velocities and a standing pattern resulting 



from the time-averaging of the forcing at large velocities. 
Analytical and numerical results will be reported featur- 
ing the simplest of these spatio-temporal behaviors. Ex- 
periments conducted with the CDIMA reaction are also 
provided which fully confirm our theoretical predictions. 
Beyond the particular chemical context that motivates 
our study, such solutions, and the conditions of their ap- 
pearance, are sufficiently generic to be applicable to a 
rather general class of pattern forming systems including 
for instance Rayleigh-Benard convection. In a sense, our 
study may be viewed as a development of the work of 
Coullet Il6| on commensurate-incommensurate transi- 
tions in nonequilibrium (spatially) forced systems. 

The experimental system under study is modelled 
within the Lengyel-Epstein scheme [fJLonce modified 
to include the effect of illumination, as |l5| : 
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Here u and v are the dimensionless concentrations of two 
of the chemical species; a, c, a, and d denote dimension- 
less parameters of the chemical system. The effect of 
external illumination is introduced through the 0-tcrms. 
This contribution can be decomposed into the mean value 
(j) and a modulation part: cj)(x,t) — 4>Q + e cos(kfX+Lut). 
For purely homogeneous illumination, e = 0, the equa- 
tions admit a solution which in the following will be 
referred to as base state: u — (a — 5<fo)/(5c), vq = 
a(l + u§)/(5uo). 

All our numerical results are obtained through inte- 
gration of the model reaction-diffusion equations 
with periodic boundary conditions by means of a pseudo- 
spectral method with a linear-implicit time-stepping. 
From here on we fix the parameters to the following: 
a = 16. c = 0.6, d = 1.07, and a — 301. These values 
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were chosen to reproduce the experimental conditions re- 
ferred to below. The large a value, corresponding to a 
strong diffusion contrast between the two species, guar- 
antees that we are far from the oscillatory regime of this 
chemical system. Consequently, only the Turing bifurca- 
tion will play a role. The remaining parameters concern 
the forcing term. One of them, (/>o, will serve as the pa- 
rameter to locate the position of the Turing bifurcation 
in the homogeneous problem. For the given parameters it 
occurs at 4>o ~ 2.3 (the base state being unstable to Tur- 
ing patterns below <fio « 2.3) and the critical wavenumber 
is k c « 1.07. 

We begin our analysis with the case of exact 1:1 spatial 
resonance k c = kf. Choosing the length L of the periodic 
domain to be 10A, where A = 2it/k c is the critical wave- 
length at instability, we fix L — 58.72 and the amplitude 
e = 0.1. 

As is known from the analysis of time-independent 
forcing, the spatial 1:1 resonance yields an imperfect bi- 
furcation to Turing patterns Q. Accordingly, the base 
state ceases to be a stationary solution and is modified 
into a non-homogeneous state for every value of the bifur- 
cation parameter 4>q. For <f>o in the stable region (above 
2.3) this state is a travelling wave (TW). The TW locks 
to the forcing wave, adopting the same wave-number and 
frequency, with only a constant phase difference. This is 
the trivial state to be expected for slow driving. For large 
uj the amplitude of the TW approaches 0. The TW's ex- 
ist to the right of the solid curve in fig. ^ 

Crossing the solid curve in fig. ^ the TW state un- 
dergoes a first instability into a state with temporally 
modulated amplitude (MTW, triangles in fig. 0. This 
is the signature of a Hopf bifurcation which introduces a 
new frequency fjj (see fig.|2t)- Note that the modulation 
of the amplitude occurs uniformly in the entire system. 
We found that for large 10 the Hopf frequency converges 
to the frequency u> of the forcing wave. 

The solutions and the transition described above can 
be rationalized in terms of an amplitude equation. Fol- 
lowing standard envelope techniques near threshold 0, 
0] , and with the forcing being invariant under the trans- 
formation t — > t+T, x — > x—uiT/kf, for sufficiently small 
e and u>, the slowly varying modulations of the travelling 
mode, in the case of perfect 1:1 resonance (kf — k c ) will 
be given by the amplitude equation 

A = fj,A - \A\ 2 A + ee?q>(-iujt) +d 2 A/dx 2 . (3) 

Using polar coordinates, A = RexpiO, we look for 
homogeneous solutions with 9 = 9o — wt. As in the 
steady case for e ^ there is a non-zero solution for every 
/i, the dimensionless distance to threshold. Its amplitude 
approaches with increasing to in accordance with our 
observations for the Lengyel-Epstcin model. 8o is the 
phase shift between the forcing wave and the resulting 
pattern. 



We further determine the stability of this solution with 
respect to homogeneous perturbations. Directly from the 
amplitude equation: 

R = fiR-R 3 + ecos(^ + 6 ), (4) 
ip = — 1 sin(^ + 9o) +w, (5) 

where we have defined: ip = + tut — Qq- Linearization 
about the locked solution (tfj = 0, R = Q = const) yields 
the following eigenvalues: A = (j, — 2Q 2 ± ^/{Q 4 — uj 2 ). 
Corresponding to fig. ^ for large u), the marginal curve 
approaches the line y, = 0, and the imaginary part of the 
eigenvalue, converges to the driving frequency, as also 
observed numerically above fli| . 

We now address the more generic case of inexact 1:1 
resonance introducing a slight wavenumber misfit, 
kt 7^ k c - To allow for a continuous variation of the mis- 
fit in a finite system, we will fix the integer wavenum- 
ber (n — kL/2ir) of the forcing to n = 10 and change 
smoothly the length L of the periodic domain. For ex- 
ample, for L = 65 the 10th wave number corresponds 
to kf w 0.97. Fig. depicts the complex behavior that 
was found changing kf from 0.9 to 1.26 (L — 70, ... , 50). 
The average illumination <fio was fixed to 2.25 during all 
of the simulations. For the purely homogeneous forcing 
(i.e. 4>' = 0) this value corresponds to a slightly unstable 
base state. 

In fig. |3 the stability domain in uj and kf is given for 
four different states. The TW state, the solution locked 
to the forcing, is the only stable solution for approxi- 
mately uj < 1. Above a roughly horizontal transition 
curve close to u> = 1.1 the TW states are unstable. For 
kf = k c and close to this point the transition occurs at 
uj « 1.15 (cf. fig.UJ and results in the MTW. 

The domain of stability of the MTW is given by the 
solid curve in fig-El Outside of this domain the attractors 
are either the TW state or one of two different new lo- 
calized states which we call soliton waves (SOW) 20] . A 
typical space-time plot is shown in fig. [2t>. Spatial plots 
of these states show that they resemble the soliton-likc 
solutions for the case of non-travelling forcing |7J but they 
exhibit the following new properties: the soliton, i.e., the 
localized suppression of the amplitude moves relatively to 
the underlying travelling pattern with a relative velocity 
which may be positive or negative. Furthermore, as the 
soliton travels along the pattern the phase of the back- 
ground pattern moves in either direction. In our case the 
integer wavenumber of the background pattern is cither 
ii = 9 or ii = 11 and the corresponding states are repre- 
sented in fig. |21 with circles and diamonds, respectively. 

Beyond the range of forcing wavenumbers shown in 
fig.|3there exist further soliton states in accordance with 
the prediction for non-travelling forcing ■ Here we will 
describe only the two states that are adjacent to the 
MTW state in fig. The states represented by dia- 
monds (resp. circles) carry a soliton that moves to the 



3 



right (resp. left). The approximate range of stability for 
these states is given by the dashed curve (resp. dotted 
curve). 

Note that the stability domains of the soliton states 
and the MTW state overlap. Remarkably, at sufficiently 
large lj the soliton states are even stable for kf = k c , i.e., 
exact 1:1 resonance. This contrasts to the case of static 
forcing Q where only the presence of a misfit can lead to 
solitons. 

The above localized states (SOW) can be described 
as quasi-periodic in the sense that their behavior can be 
reproduced as superposition of two modes, namely, the 
forcing (travelling) mode with wavelength L/10 and an 
adjacent mode with wavelength L/9 or L/ll, typically 
the closer one to the characteristic wavelength of the Tur- 
ing instability. The excited Turing mode is standing so 
that the superposition of both yields a localized envelope 
(soliton) moving to the right (L/ll) or to the left (L/9). 
Furthermore, one finds the velocity of the soliton to be 
proportional to the quotient of to and the difference of 
the wavenumbers of the two modes, which agrees with 
numerical results of the full problem. 

Note that the phase diagram depicted in fig. ex- 
hibits large regions of multistability, with the correspond- 
ing hysteresis. It is worth remarking that the stabil- 
ity boundaries here reported may well be affected by to 
finite-size effects, since possible long wave-length, phase 
instabilities may be prevented by the finite size of our 
simulations. A full envelope and phase-diffusion descrip- 
tion of this scenario in an infinite system deserves a de- 
tailed study and will be addressed elsewhere. 

A link between the SOW's and the solitons for steady 
forcing is apparent for smaller e. We decreased the value 
of e to 0.003 and determined the domain of stability for 
the soliton 11 state for u> values close to (fig-EJ. It ap- 
pears that even for vanishing u> the soliton states persist 
and the domain of stability forms a "tongue" with a finite 
range at cj = 0. For larger e and u — the SOW's do not 
exist since the driving term dominates the equations and 
forces a homogeneous pattern. Therefore it seems that 
the relative suppression of the forcing term by a effec- 
tive time-averaging due to a fast travelling forcing may 
account for the existence of soliton states for larger e. 

To prove the feasibility of controlled spatio-temporal 
forcing and to check the validity and robustness of our 
predictions in a real system undergoing a Turing insta- 
bility, we have carried out experiments on the (photosen- 
sitive) CDIMA chemical reaction, using the projection of 
travelling light patterns as the controlled forcing. Exper- 
iments were performed in a thermostatized continuously 
fed unstirred one-feeding-chamber reactor at 4 ± 0.5 C. 
Structures appear in an agarose gel layer (2% agarose, 
thickness 0.3 mm, diameter 20 mm). The gel layer was 
separated from the feeding chamber by an Anapore mem- 
brane (Whatman, pore size 0.2 mm) and a nitrocellulose 
membrane (Schleicher & Schnell, pore size 0.45 mm). 



Under the chosen set of reagent concentrations (see cap- 
tion of fig. 5), the system at dark spontaneously yields 
disordered stripe patterns with a wavelength of 0.54 mm. 
In a typical experiment, parallel light stripes with a char- 
acteristic wavelength of precisely 0.54 mm were focused 
on the gel layer and were moved in the horizontal direc- 
tion with well-controlled and constant velocity v. For 
very small values, trivial travelling stripes, following adi- 
abatically the imposed pattern, were found, as expected. 
As the passing velocity was increased they readily trans- 
formed into the modulated striped mode, as predicted by 
the theory (see fig- EK) ■ Furthermore, when considering 
a slight misfit (see caption of fig. 5), a localized structure 
was observed, propagating in the opposite direction to 
that of the stripes (see fig. Eh) -This again constitutes a 
neat confirmation of our theoretical findings. Further ex- 
periments are presently being conducted to look for other 
modes of dynamical responses to such a spatio- temporal 
modulation and will be published elsewhere. 

We have described a generic mechanism to induce new 
pattern dynamics through spatial resonance. The phe- 
nomenon is claimed to be generic for systems undergo- 
ing a Turing instability. For instance we have already 
observed the same type of response in the Id Swift- 
Hohcnberg equation. Within the framework of chemical 
Turing patterns, the consideration of the 2d case is read- 
ily accessible and points out to even richer phenomena. 
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FIG. 1: Type of attractor for different values of uo and average 
illumination <f>Q\ boxes correspond to TW solutions, triangles 
to MTW solutions (see text). The vertical line shows the 
position of the instability for homogeneous forcing. 



FIG. 2: Space-time plots of the w-component of the modu- 
lated travelling wave solution for u> = 0.5 and <f>o = 2.11 (a); 
and for the soliton solution for u) = 5, <j>o = 2.25 and kj — 0.9 
(b). 



FIG. 3: Phase diagram in the space of u) (vertical) and kj 
(horizontal), filled box - TW, filled triangle - MTW (see text 
for notation), diamond - soliton state with wave length L/ll, 
circle - soliton state with wave length L/9. We have chosen 
4>o = 2.25. The value kf = k c (L — 58.72) corresponds to the 
dashed vertical line showing the position of the perfect 1:1 
resonance. 
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FIG. 4: 'Tongue' of soliton states as a function of ui for small 
e. The dashed line approximates the boundary of the soliton 
'tongue' with n = 11 as it approaches the uj = line (spatial 
forcing of Ref.Q). Here we have used e = 0.003. 



FIG. 5: Experimental space-time plots for the modulated 
travelling wave solution (a) and for the soliton solution (b). 
The dashed line in (b) is a guide to the eye. The input con- 
centrations of reagents are 0.45 mM 12, 0.078 mM C102, 10 
mM H2S04, 1.2 mM malonic acid and 10 mM polyvinil alco- 
hol with a residence time in the reactor of 250 s. With these 
parameters the system spontaneously yields stripe patterns 
with a wavelength of 0.54 mm. (Experimental parameters: v 
= 0.13 mm/h for both cases, the imposed wavelength is 1.1 
times the spontaneous one in the case of the soliton solution 
case) . 
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